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Abstract. We analyze stationary slicings of the Schwarzschild spacetime defined 
by members of the Bona-Masso family of slicing conditions. Our main focus is on 
the influence of a non-vanishing offset to the trace of the extrinsic curvature, 
which forbids the existence of standard Cauchy foliations but at the same 
time allows gauge choices that are adapted to include null infinity in the 

evolution. These hyperboloidal slicings are especially interesting for observing 
outgoing gravitational waves. We show that the standard 1+log slicing condition 
admits no overall regular hyperboloidal slicing, but by appropriately combining 
with harmonic slicing, we construct a gauge condition that leads to a strongly 
singularity-avoiding hyperboloidal foliation that connects the black hole to ^ . 
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1. Introduction 

The principal motivation driving the development of numerical relativity has been to 
treat astrophysically relevant situations in full general relativity and make predictions 
for physical measurements — in particular for future observations of gravitational 
waves (GWs). For many years chronic instabilities in numerical codes for simulating 
the coalescence of black-hole binaries posed a major obstacle for the field, and 
correspondingly the codes currently used for black-hole simulations originate from 
a time when the main focus of the field was on the stability of evolution codes, 
rather than the accuracy of the GW content of the numerically generated spacetimcs. 
Subsequent developments have lead to the production of waveforms that are accurate 
enough for many applications in GW astronomy in the near future [1] , but succeeding 
generations of detectors (e.g. LISA [21 [3] or the Einstein Telescope [4]) will however 
significantly increase the accuracy requirements for numerical simulations. 

Following the initial breakthrough of 2005 [3 El [7] , a wealth of information on the 
coalescence of black holes has already been learned from numerical simulations (see 
[HE!] for overviews), and the accuracy of numerical codes has increased dramatically, 
see e.g. [TJ. If we are to extract the maximum possible physical and astrophysical 
information from GW observations, however, we will eventually require yet more 
accurate, and unambiguous, calculations of the GW signals from these systems. The 
calculation of the GW signal is complicated in general relativity because observers 
of signals at astronomical distances from the sources are appropriately idealized by 
quantities defined at null infinity pjjl [Til US]: i- e -> locations arbitrarily far from 
the source along lightlike directions. In most numerical codes one calculates an 
approximation to the GW signal at some set of finite distances (on the order of 
~ 100M, where M is the total mass of the system) from the source, and extrapolates 
the result to null infinity the same way one extrapolates the results from runs at 
different spacings of the numerical grid to vanishing grid spacing. This procedure 
is cumbersome, computationally expensive, and, often more importantly, prone to 
numerical errors and oddities of the numerical coordinates that dominate the error 
estimate of the final signal. 

A preferable approach would be to include null infinity in the numerical code, 
and there has indeed been much recent progress in this direction, for example in 
characteristic evolution |T3l [T4J [15] , and the hyperboloidal initial- value problem [TBI 
I17j . In this paper we are concerned with the latter approach, namely constructing 
asymptotically null slices of a black-hole spacetime, with a view to adapting our 
methods in the future to the construction of hyperboloidal multiple-black-hole initial 
data for use in numerical simulations. The construction of hyperboloidal initial data 
is however related to the method that will be used to evolve them, and in particular 
to the gauge conditions, which will be the focus of this article. 

The idea of evolving the Einstein equations as a Cauchy problem along a foliation 
of spacelike surfaces that reach null infinity was pioneered by Friedrich [HI [19] ; see 
[201 EH HH H3] f° r early overviews on numerical work. It was recognized early on that 
including null infinity in a numerical simulation naturally involves a compactification 
of an infinite physical domain onto a finite grid, and is prone to lead to singularities in 
the differential equations that we wish to solve (for a suggestion to use more generalized 
slices see |24|). Finding a convenient regularization of these equations is one of the 
major obstacles to the use of null infinity in numerical simulations. A point of view 
that has been advocated in recent years [TBI [25, 26, 27. is that before formulating the 
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concrete system of equations to be solved, the geometric structure of null infinity, which 
expresses the asymptotically Minkowskian nature of radiation spacetimes, should be 
made as manifest as possible in the coordinate gauge conditions that will be used. Such 
a strategy has worked very well for characteristic evolution [T5], [25] ■ The hope is that 
this will simplify the rcgularization problem, keep the resulting equations as simple 
as possible, and carry over many techniques from "conventional" numerical-relativity 
simulations on finite Cauchy slices. 

Before doing this, one must settle on gauge conditions. It is common practice 
in numerical relativity, and ubiquitous in current black- hole-binary codes, to evolve 
gauge conditions with hyperbolic equations. For example, the popular moving- 
puncture method employs gauge conditions chosen from the Bona-Masso family [29] 
(see however (30] [HUGS] f° r elliptic- hyperbolic approaches). 

Although one ultimately wants to simulate a dynamical spacetime, stationary 
solutions play an important role: in many applications perturbations ultimately 
radiate away, and a stationary solution is approached. This is true in particular for 
astrophysical black hole spacetimes. It is therefore fruitful to first think in terms of 
stationary solutions with respect to the gauge conditions. In fact, it is somewhat 
surprising that the search for stationary representations of black holes which arc 
compatible with the gauge conditions one is using was not posed much earlier as one of 
the fundamental questions in numerical relativity. Stationary solutions also motivate 
more suitable numerical methods and forms of initial data. The importance of time- 
independent descriptions of Minkowski-spacetime as a starting point for numerical 
approaches to the hyperboloidal initial value problem and as a test-case for the 
conformal field equations approach has been discussed in [23] . For recent work on the 
associated problem of freezing the coordinate postion of null infinity in compactificd 
evolutions see [33l [34] . 

In [35l [36[ [37] we have discussed explicit stationary representations describing 
a nonspinning black hole, which are consistent with the specific gauge conditions 
used in the moving-puncture method and represent a trumpet geometry: the slice 
extends from a throat at some finite value i?o of the Schwarzschild radial coordinate 
out to spatial infinity. In [31] such maximally sliced "trumpet data" have been 
matched explicitly to a family of stationary hyperboloidal slices. In this paper, 
we construct such "hyperboloidal trumpet" slices directly from an evolved gauge 
condition, with the hope that it will provide a useful starting point for developing 
a method for the simulation of black-hole-binary spacetimes with asymptotically null 
slices. We first study the compatibility of the Bona-Masso family of slicing conditions 
|29j with regular stationary hyperboloidal slices of Schwarzschild-Kruskal spacetime. 
We find incompatibility for the popular "1+log" subfamily, but compatibility for 
harmonic slicing and certain "hybrid" conditions, which we construct to yield 
stationary slices which are both hyperboloidal and singularity-avoiding, and which 
are therefore interesting regarding the generalization of the moving-puncture method 
to the hyperboloidal problem. 

2. Preliminaries 

2.1. Geometry of hyperboloidal slices 

Before embarking on technicalities, let us note that a crucial geometric quantity in 
our work is the (trace of the) extrinsic curvature — in this paper we will choose 



Stationary hyperboloidal slicings with evolved gauge conditions 



4 



the convention that the extrinsic curvature of a 3-dimcnsional spacclikc hypersurfacc 
is defined as the Lie derivative with respect to the timclikc unit normal n a of the 
induced metric h a b, = +hC n h a i,, with a positive sign. This sign convention is 
unfortunately not common in numerical relativity, but we believe it is beneficial for 
physical intuition: positive mean curvature K = h ab K a b then signifies an expanding 
volume element, while negative K signifies contraction. 

While we arc interested in the asymptotic structure of slices, in order to simplify 
notation we will actually not directly work with the conformal compactification picture 
[TU] . which has become a standard tool [35], but we nevertheless want to define 
hyperboloidal slices as regular spacelike hypersurfaces in a compactified spacetime in 
the sense of Penrose [10] . It then follows for the physical trace of the extrinsic curvature 
K, that K > at future null infinity J? + , and K < at past null infinity Jp^ . The 
name "hyperboloidal" stems from the fact that such surfaces are analogous to the 
standard hypcrboloids t 2 — x 2 — y 2 — z 2 = (3/K) 2 in Minkowski space, which provide 
the simplest example. A wider class of important standard examples are provided 
by the constant mean curvature (CMC) slicings, defined as K — const., which are 
important for cosmology and very well studied in spherical symmetry [39] [40] [41] . For 
the purpose of comparing our results to CMC slicing, it is sufficient to recall that the 
stationary lapse a for a constant K = K in spherical symmetry reads 

I fKR 2M 

where R is the areal radius, M the mass of the Schwarzschild black hole and C an 
integration constant that is of no interest for us. 

The key difference between a Cauchy slice and a hyperboloidal slice of an 
asymptotically flat spacetime, for the purpose of this paper, is in the asymptotic 
behavior of the standard lapse function a and shift vector [3 a [42] . It is well known 
that for a Cauchy foliation one has that a — * 1 and \(3 a \ — > as r — > oo, while for 
hyperboloidal foliations we have that 

a = (D{r), \p a \=(D(r 2 ) for r -> oo, (2) 

where r is a radial coordinate that asymptotically behaves like the Schwarzschild 
radial coordinate. For a stationary hyperboloidal slice that reaches out to J? + , where 
the mean extrinsic curvature asymptotes to a positive value, the shift vector points 
inward at large separation, while at J?~ it points outward. This is intuitively clear 
from looking at a Penrose diagram such as figure [TTt in order to keep say, at 
a constant coordinate location, the shift vector has to point along J^ + , i.e., along 
an ingoing null surface. Likewise, to keep a black hole horizon in place, an outward 
pointing shift vector is required, and consequently the shift vector must be expected 
to change its sign for stationary slices that connect a black hole to (or a white hole 
to ^~). Some attention to the direction of the shift vector will be required below. 

2.2. Bona-Masso slicing conditions 

In 1994, Bona et al [53] proposed a family of slicing conditions that are formulated 
as evolution equations for the lapse a. Introducing a positive, otherwise arbitrary 
function f(a), these conditions can be written as 

{d t - C )a = a 2 f{a)(K - K ) , (3) 
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where Cp denotes the Lie derivative along the shift vector (3 a . For / = 1 we recover 
harmonic slicing, / = n/a, n £ K is called 1+log slicing (with n = 2 the most popular 
choice for black hole evolutions); maximal slicing can be viewed as corresponding to 
the limit / — * oo (see e.g. our work [37]). The general properties of these slicing 
conditions have been studied in much detail [23 HH HH US [^1 [47l [37l [35] and they 
have been implemented in numerical simulations with great success, in particular in the 
breakthrough simulations that have established the moving puncture method [6[ [7] . 
So far these conditions have been used in the context of asymptotically Euclidean 
foliations, and the function Kq has usually been taken as a constant, and explicitly or 
implicitly been set to zero. 

In this paper, we ask the question of whether the slicing condition (J3]) is also useful 
in the hyperboloidal context. We will still consider Kq a constant for simplicity, but we 
will now take it seriously. We find that the value of Kq indeed has a strong influence 
on the stationary foliations of the Schwarzschild spacetime, which for Ko ^ differ 
considerably from the stationary solutions found in [?7J [37J [351 [35] . 

Note first that for Kq = stationary solutions will satisfy 

Cpa = -a 2 f(a)K , (4) 

which is indeed consistent with a — > 1, (3 a — > 0, K — > at infinity, i.e., the condition 
is consistent with asymptotically Euclidean slices. Looking at the right hand side 
(RHS) of ([3]), one finds that Kq > locally has a decreasing effect on the lapse, while 
Ko < has the effect of increasing the lapse function. The question arises whether 
Kq 7^ is compatible with stationary slices of the Cauchy type. At first sight one 
might naively guess that Kq > leads to asymptotically constant-mean-curvature 
(CMC) slicings that reach J^ + , whereas Kq < leads to asymptotically CMC slicings 
that reach J^ - . We will see that this is indeed the case, when regular CMC slices 
exist. For an asymptotically Euclidean slice, however, the RHS of ([3J approaches a 
constant if Kq ^ 0, while the left-hand side must go to zero, and therefore a stationary 
asymptotically Euclidean slice is incompatible with the choice of Kq ^ 0. We will now 
investigate the construction of asymptotically CMC slices. 



3. Construction of stationary hyperboloidal solutions 

3.1. 1+log Slicing 

We start our detailed analysis with the general version of 1+log slicing and n = 2, 

(fit - Lp)a = 2a(K - K Q ) . (5) 
After introducing our notation, we shall show explicitly that for Kq there are no 
overall regular time-independent solutions of ([5|) in spherical symmetry that lead to a 
Cauchy foliation or become asymptotically CMC slices. 

The calculations are restricted to the spherically symmetric spacetime of a 
Schwarzschild black hole with total mass M, 

ds 2 = - (l - dT 2 1 dR 2 + R 2 dCl 2 , (6) 

where T and R denote the standard Schwarzschild coordinates and dQ = dip 2 + 
sin 2 8 d9 2 . In the case of stationary (i.e., time-independent) foliations, the trace of 
the extrinsic curvature can be expressed as |39j 
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where ' denotes the derivative with respect to the areal radius R and 



(3 = f3 R /a = ±y/im ■ (8) 
Furthermore, the lapse and the shift are related through 

a?-f = l- 2 -f. (9) 
Thus, time-independent solutions of §5§ satisfy 

4/3 

Pat = -2{K - K ) = 2/3' + -£ + 2K (10) 



R 



4i?a 2 + 6M + 2R ( K Ra/ a 2 + 2M - 1 - 2 

■ (11) 

R (Ra 2 - 2Ra + 2M-R) K ' 

Note that (3 was defined as the positive root of © to obtain (pTTj) . We could as well 
have chosen the negative sign, which has the same effect in (jlll) as changing the sign 
of Kq. For clarity, we make the following statements for the class of solutions that 
have positive shift for small R. However, it should be kept in mind that swapping the 
sign both of (3 and Kq has no effect on the solution for the lapse. In |Appcndix A| we 
also write down a coupled system of differential equations for lapse and shift, which 
is equivalent to ©, and where an explicit tracking of the sign of the shift is not 
necessary. 

By closely investigating (|10[) we immediately convince ourselves that the 
stationary lapse can neither go to unity, nor can it be CMC-like in the range of 
large R. For a — ► 1, we find that a', (3 and (3 1 vanish in the limit R — *■ oo, which 
contradicts p0[) for all K 7^ 0. If a ~ R and a' — * const., as in the case of the 
stationary CMC solution ([1]), we can conclude [3a' ~ R, whereas the RHS of (fTU]) 
approaches a constant. The only power law that does not lead to a contradiction is 

a~\/# => K~1/Vr (12) 

and below we will indeed find a family of solutions that show this asymptotic behavior 
for large values of R and positive Kq. 

However, let us first discuss solutions with negative Kq (and positive shift). We 
analyze the RHS of (|11| by noting that the denominator vanishes along the curve 

R _ 2 
M ~ 1 + 2a - a 2 ' 

that is characterized by the limit a — > 1 + y/2 for R — > 00. The numerator of (jTTJ) 
vanishes for 



(13) 



3M K 2 R 2 K Q y/R 
~2R + ~ 8 



K 2 R 3 + 8M, (14) 



which grows unboundedly for R — > 00 and Kq < 0. For every point (a, R), we can 
easily determine the sign of a' by investigating in which region it is located between 
the functions (fl"3|) and (fT4"|) . If a curve that solves the differential equation (jTTJ) enters 
a region with negative slope, it is already clear from figure [1] that it will be driven 
towards the dashed curve (|13|) . either directly in region (IV) or by going from (I) to 
(II). Note that from the discussion above we know that a' cannot vanish in the limit 
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Figure 1. The analysis of stationary states that satisfy the 1+log equation ((SJl 
with negative Kq (here, Kq = —0.5/M). a' > is indicated by regions (II) and 
(III), the dark gray regions (I) and (IV) illustrate a' < 0. Their boundaries are 
given by (1 1 3 I t (dashed curve), H14I I and a = \/l — 2M/R. All solutions become 
singular where the dashed curve is hit. Some examples which we have calculated 
explicitly are shown for clarity. 



R — > oo. One might think that solutions can "escape" to infinity in region (III), but 
for sufficiently large a and R, equation (|lip reveals 



so that for further increasing R and a, all solutions coming from (III) must enter (IV) 
by crossing the line (fTl)) . 

This finally leads us to the conclusion that every integral curve which satisfies 
pip will inevitably hit the function (|13|) and become singular. Put another way, for 
negative Kq there exist no smooth solutions of (fTT|) that are well defined for arbitrarily 
large R. For completeness, we note that the picture changes for Kq < — 2.5896/M, 
since then the functions (fT5|) and (jT5J) do not intersect at positive R and a. However, 
our arguments are based on the limit R — > oo and therefore still hold for this case. 

Let us proceed to positive Kq. It is easy to prove that now, the curve of vanishing 
numerator (|14p approaches the limit a — > 1 for R — > oo. Thus, there exist solutions 
of (fTTjl that do not become singular but are well defined for arbitrarily large R (and 
a > 1 + v2)- We confirmed numerically that these are the predicted solutions with 
a ~ 2\[Kq~R (R — > oo). Since /3 does not change its sign, the slices are characterized 
by a negative K and thus connect the black hole to J^~. A visualization of our results 
is given by figure [21 One can also observe that several integral curves which cannot 
grow unboundedly, hit the boundary of the domain, a = yl — 2M/R. At this point, 
we find (3 = 0, so that (JTUJ) is not a well-formulated differential equation any more. 
However, by exploiting the relation ([5]), we find 




(15) 




(16) 
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Figure 2. The equivalent of figure [T] but with positive Ko (the chosen example 
is Kq = Q.l/M). Note that for increasing R, the solutions of JTTJ are driven away 
from the dashed curve, to eventually escape to infinity or change the sign of /3 (see 
text). The included log- log plot illustrates how the outer left curve approaches 
a = (straight dashed line) for large R. 



for a — > y/l — 2M/R. Hence, the shift changes its sign and becomes negative. The 
curves that hit the boundary of the domain can therefore be smoothly extended by 
allowing (3 to cross zero at this point. For negative shift, (fTTj) has to be changed as if 
Ko i ► — Ko. We therefore conclude that all these solutions run into a singularity at 
finite R. Interestingly, the unique solution that is defined as the one that crosses the 
common root of (fT5|) and (|14[) (see for example [37J EHl EH] for the discussion of the 
standard 1+log case) is of this type. Figure [3] shows this particular example, which 
we refer to as the "critical solution" . As we see, due to the discussed properties of 
(llip there is no overall regular critical stationary solution of 1+log slicing with offset 
Ko^O. 

3.2. Harmonic Slicing 

In this section, we study the harmonic slicing condition 

(d t -Cp)a = a 2 (K-K ) . (17) 

We find that in this case, K ^ does lead asymptotically to CMC solutions at large 
distances from the black hole. As expected, care has to be taken in assuming the sign 
of the shift, since for positive Ko it changes from '+' inside the black hole to '— ' in 
the limit R — *■ oo. 

The analysis of harmonic slicing can be motivated by considering the generalized 
Bona-Masso slicing condition (J3j) , restricted to f(a) = a K . First, k may be any 
exponent. Looking for stationary solutions that may asymptotically approach the 
CMC solution |T|), we exploit the assumption a ~ |A'|i?/3 (R — > oo) to determine k 
and K. Starting from ([3]), we find 

pa' =-a K+1 (K - Ko) (18) 
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Figure 3. The lapse, the shift and the trace of the extrinsic curvature of the 
numerically-calculated critical stationary solution of 1+log slicing with offset 
Kq = 0.1/M. For small R we recover a behavior similar to the stationary trumpet 
solution of standard 1+log 1371 . At R fa 3.76M, however, f3 changes its sign and 
K seems to settle close to Ko, before finally all quantities become singular at 
R « 4SM. 



\K\KR , (\K\R> 



9 -\Kq-k) y-^-j ■ d9) 

By comparing the exponent of R on both sides of the estimate, we conclude that k = 
and 

K^'^ (fl-oo). (20) 

Vanishing n obviously leads to the slicing condition (jTTJ) . 

We now discuss time-independent solutions of (JTTJ) in analogy with our analysis 
of the 1+log condition in the previous section. We note that stationary states with 
positive shift satisfy equivalently to (| 1T[) the following ordinary differential equation, 



a 



3M + R{2a 2 - 2 + K RJa 2 + ±f- - 1) 



(21) 

(2M - R)R y ' 

The set of points of vanishing numerator for R > 0, a > is still described by 
(|14| . The denominator is zero at the horizon, R = 2M. In this special case, we 
can analytically express the regular singular point (i? s ,a s ), at which both functions 
intersect, as 



a s = i y 1 + 2K 2 - 2K ^Kl + 1 , i? s = 2M . (22) 

Note that we face a different situation than in the 1+log case. Here, the curve of 
vanishing denominator is a vertical line in the a-R diagram. Therefore, every solution 
of (f2Tj) that smoothly connects the inside and the outside of the black hole must 
cross the line of vanishing denominator. To prevent the integral curve to become 
singular there, the numerator should vanish at the same point, i.e., the unique regular 
stationary lapse must pass the point given by (|22p . 
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Figure 4. The stationary lapse (left) and shift (right) of harmonic slicing 11171 , 
Each curve is obtained by the unique solution of l|21j l that passes the regular 
singular point II22H , The following cases are displayed (from top to bottom) 
MKq = —10, —I, —0.1, —0.01, 0, where vanishing Kq is shown as a dashed line. 




01 0.1 1 10 100 1000 10 4 

Figure 5. The (negative) trace of the extrinsic curvature of the stationary states 
shown in figure [4] The horizontal dashed lines correspond to 3Ko/2, the long- 
dashed curve illustrates the analytically known case of vanishing Kq . 



We analyze the resulting critical solution for the lapse by numerically integrating 
(|2Tj) . Our results show that K indeed asymptotically approaches 3K /2 for large 
R in the case of both negative and positive Kq. The plots for different choices of 
Kq are shown in figures [5] and [7j Furthermore, the left panels of figures [4] and [6] 
confirm that the lapse asymptotically grows unboundedly similar to CMC slicing, i.e., 
a~ \K\R/3~ \K \R/2. 

For negative Kq, the shift f3 is decreasing for small R until the A'o-dcpcndcnt 
minimum is reached, then the slope changes its sign and (5 — > oo for R — > oo. Thus, 
the shift is positive everywhere. The trace of the extrinsic curvature remains negative 
for all R. We conclude that the stationary foliation of harmonic slicing with negative 
offset connects the singularity R = (at infinite proper distance) to past lightlike 
infinity <?~ . 

Positive offset, on the other hand, always leads to stationary states where a hits 
the boundary of the domain, a = y/T— 2M/ R. (3 vanishes at this points and changes 
its sign from positive to negative, as described for 1+log in the previous section. By 
similar calculations one easily proves that again 



= -Kq . (23) 



Q =v 1 — W 

The numerical calculations also suggest that the root of the shift converges 
monotonically to the horizon for increasing Kq, which may lead to problems when 
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Figure 6. The stationary lapse (left) and shift (right) of harmonic slicing l|17|l 
with non-negative offset MKq = 1.25,0.25,0.05,0.01,0 (from top to bottom on 
the left panel, vice versa on the right panel; Kq = as dashed line). 




Figure 7. The trace of the extrinsic curvature of the stationary solutions that 
are displayed in figure [6] The dashed lines correspond to 3Kq/2. 



numerically solving (|21|) for large Kq. However, as a proof of principle we show 
several choices Kq < 1.25M" 1 in figures [6] and [71 The important difference to the 
solutions with negative offset is that here, changing the sign of f3 induces a change of 
the sign of K. Asymptotically approaching a constant positive K yields a slicing that 
smoothly connects the singularity of the black hole with future null infinity J? + . 



3.3. Combined approach 

Let us summarize what we have found in the previous sections. Starting from a 
general form of Bona-Masso slicing ([3]), we analyzed the 1+log case (/ = 2/a). After 
concluding that in this case there is no overall regular stationary state that connects 
the black hole to we found that harmonic slicing (/ = 1) is suitable for having 
asymptotically CMC-like stationary states. The drawback in this case is that harmonic 
slicing is only marginally singularity-avoiding — the slicing comes arbitrarily close to 
the singularity of the black hole. What we propose now is to combine both attractive 
properties, the strong singularity avoidance of 1+log and the asymptotically constant 
K > of harmonic slicing (with positive offset) by simply adding both conditions in 
an appropriate way, 

/ = - + 1 (24) 
a 

=> (d t -Cp)a = (2a + a 2 ){K - K ) . (25) 
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Figure 8. The throat (root of the lapse) of the critical stationary solutions of 
(125 I t - The horizontal dashed lines correspond to the limits Ro — ♦ 2M for Ko — > oo 
and R -> 1.3955M (see text) for K -> 0. 



The function / is obviously dominated by the 1+log part for small a (i.e., inside the 
black hole) and by the harmonic part for large a (which is in the case of CMC slicing 
equivalent to large distances from the black hole). 

We analyze stationary solutions of (|2"5")) in the case of a Schwarzschild spacetime 
by using the established methods and equations [recall especially ([7]), ([5]) and JSJ], 
The time-independent form of (|25p can be written as 



(2 + a) MI + R{2a 2 -2 + K RJ a 2 + if- - 1 



R[2M - R(l + 2a)] ^ 

The critical solution we consider is defined by passing the common root of the 
numerator (sec (1141) that is valid for arbitrary /) and denominator 

R 2 2M-R 

M=TT^ * a = ~2R-- (27) 
Before we further investigate the set of solutions for positive Ko, we give some exact 
analytical expressions for the case Kq = 0. Under this assumption, time- independent 
solutions of ([25]) satisfy 

&ol = (2 + a) (V + |T 

- ,*!L = [(p + *: ]dR 



a J \f3 R 
2 2AI C 2 (2 + a) 2 

The (positive) integration constant C 2 is determined by picking the critical solution 
as described above. As the regular singular point, we find 

Rs = ^ + - M w 1.535M , (29) 
0.151 , (30) 



13-3 



4 

and the constant is given by 



C 2 = — ^13a/13 - 35J M 4 « 0.391M 4 . (31) 
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Figure 9. The stationary lapse (left) and shift (right) of our combined slicing 
condition ((25} with non-negative offset MK = 1.25, 0.25, 0.05, 0.01, (from top 
to bottom on left panel, vice versa on the right panel; Kq = as dashed line). 
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Figure 10. The trace of the extrinsic curvature of the stationary solutions that 
are displayed in figure [9] The dashed lines correspond to 3i\o/2. 



The implicit solution (|28]) is now uniquely determined and can also be written explicitly 
as roots R(a). Note that the roots have to be picked carefully in order to obtain the 
smooth, monotonically-increasing function we are interested in (see also for a 
similar discussion) . One important feature is that the resulting lapse starts at Rq > 
with a(i?o) = and a'(R n ) > 0. Hence, the behavior of strongly singularity-avoiding 
conditions is retained and, expressed in the terminology of [37], the region close to 
the black hole is a "trumpet" with the throat Rq « 1.3955A/. However, the case 
Kq = reveals a — * 1 for R — * oo, and is therefore not qualitatively different from the 
standard 1+log. 

If we now allow Kq to be positive, we find by numerically integrating (|26[) that 
a still starts at the throat Rq(Kq) (see figured]) but then follows the behavior that 
was already observed for harmonic slicing, which means that, at finite R, the first 
positive shift changes its sign and the lapse asymptotically approaches the CMC limit 
a ~ KqR/2. The trace of the extrinsic curvature changes sign as well (from negative 
to positive) and approaches the predicted limit 3A'o/2 for R — > oo. Our results are 
displayed in figures [9] and [TUJ 

In order to point out the new quality of the proposed slicing, we also calculate 
the embedding in a Carter-Penrose diagram. For this purpose, we transform our 
integrated quantities a(R) and (3(R) to standard Schwarzschild coordinates using 
the height function approach [311 [SD]- Transforming these coordinates to Kruskal 
coordinates and compactifying them is straightforward and was already described in 
[37] . The resulting diagram is shown in figure [TT1 As discussed before, the slices 
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Figure 11. The Carter-Penrose diagram of the slicing defined by the critical 
stationary solution of (1 25 I t and Ko = 1M~ 1 . The dashed line illustrates the curve 
of constant R = Rq pa 1.8831M. All slices connect the throat Rq to Each 
displayed time step is 5M. 

indeed go towards J? + (recall, positive K) and approach timelike infinity i\ along the 
throat, as it is typical for strongly singularity-avoiding slices. 

4. Discussion 

Let us summarize our results. 

For the 1+log slicing condition used in moving-puncture simulations of black-hole 
spacetimes, we have previously found that it is possible to find a stationary solution 
for a Schwarzschild black hole. This solution represents a trumpet geometry: the slice 
extends from a throat at some finite value i?o of the Schwarzschild radial coordinate 
out to spatial infinity [35, 36, 37]. In the present paper we have extended these results 
to spacelike slices that instead reach null infinity, and which are commonly referred 
to as hyperboloidal slices. Such hyperboloidal slices asymptotically approach a finite 
value of the mean extrinsic curvature, and include the constant-mean-curvature slices. 
One potential approach to produce trumpet slices that also extend to null infinity is 
to modify the standard 1+log slicing condition to include an offset term that leads 
to constant-mean-curvature at large distances. We have shown that there are no 
regular slices that satisfy this condition, given by ([5]). Also, such an offset term 
is an obstruction to construct regular stationary Cauchy slices for the 1+log slicing 
condition. 

If we instead deal with the analogous modification of harmonic slicing, 
equation (|17[) . regular stationary slices do exist. However, harmonic slices reach 
the black-hole singularity, and for moving-puncture-like simulations we would prefer 
singularity-avoiding slices. Both requirements singularity-avoiding 1+log slicing 
near the black hole and CMC harmonic slices at null infinity — can be met by using 
a hybrid slicing condition of the form suggested in (|25[) . 

(dt-C p )a= {2a + a 2 )(K- K ) . 

We have demonstrated the efficacy of this approach by constructing slices of the 
Schwarzschild spacetime that extend from Rq (1.3955M < i?o < 2M, dependent 
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on the choice of K ) to null infinity. Both for the hybrid condition and harmonic 
slicing, the sign of the trace of the extrinsic curvature (and therefore the offset K ) is 
directly related to the outer end (R — > oo) of the slices via 

„ Z T , ,. . 4 f J+ for K > 

A — > -An slices extend to 4 „, r . 

2 \ for A < 

There is obviously a large freedom in constructing gauge conditions and stationary 
slicings with similar properties, and in choosing (or generalizing to a function) the 
constant Kq. What works best will presumably depend to some degree on the 
application. 

The results presented in this paper are intended as first steps toward producing 
multiplc-black-holc "hyperboloidal trumpet" data. Numerically evolving such data 
would require a sufficiently regular treatment of null infinity in Einstein's equations. 
The philosophy we advocated in the introduction was that the most practical approach 
may be to perform the regularization of Einstein's equations after choosing the type 
of initial data to be evolved, and the gauge conditions to be employed. As such, we 
suggest the data and gauge conditions we have presented here as one possible starting 
point for deriving a method for null evolutions that is also close in spirit to current 
moving-puncture simulations. 
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Appendix A. Bona-Masso slicing stationary equations as a coupled 
system for lapse and shift 

In the main text we have used a single differential equation for the lapse and the 
quadratic algebraic equation ((9]) to construct our slicings — tracking the solution of 
([5]) with the appropriate sign of the shift vector. For completeness, we now also write 
down a coupled system of differential equations for lapse and shift, where an explicit 
tracking of the sign of the shift is not necessary. We write the metric of a Schwarzschild 
black hole as, 

ds 2 = - [a 2 - 7 (/3 fl ) 2 ) dt 2 + 2 1 (3 R dtdR + jdR 2 + R 2 dQ 2 , (A.l) 

where t is the time coordinate of the stationary slice, R the standard Schwarzschild 
radial coordinate, d$l = def) 2 + sin 2 9 d6 2 , and a, (3r and 7 are functions of R. One 
finds from the Einstein equations that 7 = 1 /a 2 . Note that we now use the radial 
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shift component j3 R instead of (3, in order to keep track of the direction of the shift. 
In the case of time-independent foliations the extrinsic curvature reads 

(r*(/? fl ) 2 7)' 

where again ' denotes the derivative with respect to the arcal radius R. The Einstein 
equations together with the stationary Bona-Masso slicing condition Q then imply 



, _ f(a) (K r 2 j3 R + 2ro? + (3 - 2r)a) 



ra(ra(a — f(a)) — r + 2) 
and 

f(a) (2K a ra 3 (2ra 2 - r + 2) + a 2 f3 R (7ra 2 - 3r + 4) + r {f3 R ) 3 ^j - 2a 3 (3 R 
^ 2ra 3 (ra(a — f(a)) — r + 2) 

where r = R/AI. 



(A.3) 
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